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Abstract 

The effect of the joint drive system with gear reduction for a generic 
two-link system is studied. It Is done by comparing the kinetic energy of 
such a system with that of a direct drive two-link system. The only 
difference are two terms involving the inertia of the motor rotor and gear 
ratio. Modifications of the equations of motion from a direct drive system 
are then developed and generalized to various cases encountered in robot 
manipulators. 

Introduction 

Formulating the equations of motion for a robot is an Important part of 
robot analysis that will provide necessary information for the design of 
control laws and mechanical components. Before the process of formulation can 
begin, idealization of the robot system into a model amenable to analysis has 
to be performed. Assumptions such as rigid bodies, perfect revolute joints, 
complete isolation between electrical phenomena and mechanical motions and 
idealized torque transmission in the gear trains are commonly made [1-5]. By 
removing one or several of these assumptions, one can come up with models with 
different levels of fidelity. The price to pay is the increased complexity in 
the resulted equations and possible numerical difficulties. But sometimes, 
the price has to be paid in order to obtain equations that correspond better 
with the important characteristics of the actual robot dynamics. In this 
paper, the effect of some idealizations of joint drive systems in the commonly 
used model will be investigated. 

Although many robot joints are driven by motors through the use of gears 
with reasonably high (on the order of hundred) reduction ratio, the commonly 
used model does not include any detail of the drive system. Strictly 
speaking, the model of the generally used multi-body system is only directly 
related to a direct drive robot. For this model to be applied to a robot with 
torque transmission and amplification, certain rules are usually implied. It 
is generally believed that the torque at the joints are equal to the product 
of the corresponding motor torques and reduction ratios. It is also known to 
some that the rotational inertias of the motor rotors when amplified by the 
corresponding reduction ratios squared should be included in the link inertia 
to resist motion. These implications can be found in textbooks on automatic 
control, such as [6], but they are rarely spelled out in robot literature and 
their validity is not established. 
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In this paper, a simple model will be used to study the effect of gear 
reduction on system dynamics. Specifically, the necessary modifications to 
the equations of motion for the commonly used simply jointed robot model will 
be presented. An example will also be used to demonstrate the effects of gear 
reduction to the equations of motion. 

System Model 

In the present study, It Is assumed that the motor rotor Is the only 
massive element In a joint drive system that will contribute to the modifica- 
tions of the equations of motion from that of a multi-body direct drive 
system. To study this problem, one might be tempted to just study a par- 
ticular system by Including the rotors In the model and deriving equations for 
the system. Simulation can then be performed to hopefully reveal the effects 
caused by the Inclusion of rotors in the model. This Investigation, however, 
will be only specific for the particular system simulated and will not shed 
too much light for a system with very different geometry and mass distribu- 
tion. The methodology adopted here Is to understand the general effects of 
the rotors and thus to enumerate the suitable modifications of the equations 
due to them. In place of a specific model, a generic system model should be 
used and the resulted equations studied to Identify and generalize the effects 
of rotors. 

The generic model chosen is shown in Fig. 1. Body C Is the carrier of 
the motor whose rotor together with the connecting shaft and attached gear 
forms a rigid body R. The driven link D has an attached gear which meshes 
with the gear In R. This is an Idealized system of two links with simple 
drive system. Since It Is assumed that the rotor is the only massive element 
In the drive system, this model Is sufficiently comprehensive for the study. 
Instead of letting C be a link jointed to the base, it Is allowed a general 
motion relative to the base. This Is an Intentional choice so that C can be 
any link of a multi-link system. The system can therefore be considered as a 
subsystem of an overall system. It can be seen that the linkage in Fig. 1 
Involves a closed-loop topology and formulations for the commonly studied 
open-loop multi-link system cannot be applied here. In particular, Newton- 
Euler formulation Is not very convenient for this system and is especially 
difficult to generalize. Lagrange's or Kane's formulation is more suitable 
for the present study because generalized active and inertia forces can be 
considered as being contributed from individual elements of the system. The 
overall system can be thought of as having n generalized coordinates, qi, ..., 
q n , and q, the joint angle between C and D, can be one of them. 

Formulation of Dynamic Equations 

Vector-dyadic formalism will be used in all the following formulations. 
Here, a vector and a dyadic are defined as abstract entities which are 
invariant with respect to unit vector bases [7-9]. They can be expressed in 
terms of any unit vector basis or bases but some operations among them can be 
reduced without being expressed in any basis. Vectors and dyadics as used in 
the following are therefore different from column matrices and square matri- 
ces, which are just congregates of numbers. However, when the vectors and 
the dyadics are represented in the same vector basis, their operations can be 
facilitated by matrix operations. 
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The first step In Lagrange's formulation Is to derive the kinetic energy 
of the system. For the subsystem in Fig. 1, the contribution to the kinetic 
energy Is 

K = i X m( N v P ) 2 d) 

c C+R+D 

where N Is an inertia reference frame, m and N v P are, respectively, the mass 
and the velocity in N of a generic particle P in the system, and 2, denotes 

C+R+D 

sumnatlon over all particles in bodies C, R and D. In the notation for a 
velocity, an acceleration, an angular velocity or an angular acceleration, a 
left superscript Is used to denote the reference frame the quantity Is referred 
to. In the sequel, when the left superscript Is omitted In any of these nota- 
tions, reference frame N Is Implied. Applying the theorem called one point 
moving on a rigid body [9], one can express the velocity of a generic particle P 
of R or D by 


P C C P 
v = v + v 


( 2 ) 


where C v P Is the velocity of P In C, and v C Is the velocity of l, which is a 
point of C that coincides with P at the instant under consideration. 
Substitution of equation (2) In equation (1) yields 


K = l f l m(v P ) 2 + l m[( C v P ) 2 + 2 v C • C v P ]} (3) 

c E R+D 

where E Is a fictitious rigid body that moves exactly like C but has exactly 
the same mass distribution as that of C, R and D altogether at the Instant 
under consideration. 


With the application of a kinematic theorem for two points fixed on a 
rigid body [9], one can express the velocity of a generic particle P of E as 


v p - 


C 

(!) X 


.OP 


(4) 


where w c Is the angular velocity of C in N, Q Is an arbitrary reference point 
on C and rQ p is the position vector from 0 to P. The use of equation (4) 
brings the first term In the right hand side of equation (3) to the form of 

Y, m(v P ) 2 = m F (v ^) 2 + w c • 1^ • w c + 2m f v^ • (» c x r^ ) (5) 

E t 


where 

I E/Q - l m [(r 0P ) 2 U - r 0P r° P ] (6) 

E 

* OE* 

and, mg and E are the mass and the mass center of E, respectively, while r 

Is the position vector from Q to E*, and II is a unit dyadic. It should be 


299 



noticed that E and I E/ Q are not fixed In C. As for the second term In 
equation (3), the following equations can be similarly derived by suitably 
choosing a reference point. Firstly, 


i m( c vV 

R 


mW ♦ V 


jR/R* . C # R 


( 7 ) 


where 


jR/R* 


X m[(r R * P ) 2 U - 
R 


r R*P r R*P-, 


( 8 ) 


and, mp and R^ are the mass and the mass center of R, while C v R and C w R are the 
velocity of R In C and the angular velocity of R in C, respectively, and r R * p 
Is the position vector from R to P. Next, 


X m( c v p ) 2 . C « D . I D/Q ‘ . C W D 
D 


(9) 


where 


jD/Q' 


X m[(r Q,p ) 2 U - 
D 


r 0'P r 0'P-| 


( 10 ) 


QP CD 

and, r and u are the position vector from Q' , a point on the joint axis 
to P and the angular velocity of D In C, respectively. Then, making use of’ 
equation (4) for v p , one can write 


X m v C • C v P = v^ • m D C v R + W C ♦ c h R/Q 
R R 

where 

W'O - l m r QP x V 
R 

C R/Q 

and H is the angular momentum of R about Q in C. Similarly, 

X m v C • C v P = v^ • m n C v D + w c • 

D D 

where 

C H 0/Q - z m r« p X c v p 
D 

Since R* Is fixed in C which implies 


( 11 ) 


( 12 ) 


(13) 


(14) 


C , R 


* 


= 0 


(15) 
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and 


r r 

L w K = y q a 

where y is the gear ratio, it follows that 


c H R/() .Mia 


and 


C. R T R/R* . C R 2 , -2 


V' • I 


w 1 ' = \f J q‘ 


where J is the axial moment of Inertia of R. Also since 


C D 

w = q c 


equation (9) can be rewritten as 


l m( C v P ) 2 = q 2 (c • I D/QI • c) 
D 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


Substitution of equations (4) - ( 18) in equation (3) yields 


K = \ {m E (v Q ) 2 + w c • I E/Q • w c + q 2 (p 2 J + c • I D/Q ' • c)} 


+ w c • (r QE * x m E v Q + p J q a + V, H L " v ) + • ’V 


C..D/Q. 


0 CD* 

1 1* • u 


( 21 ) 


It can be seen In equation (21) that the kinetic energy of the system Involves 

only two terms, \ p 2 Jq 2 and yJq w c • a that depend exclusively on the motor 
rotor R. Other contributions to K due to R are lumped in those terms 
involving the fictitious body E. 

Consider now another system S' consisting of only two links C 1 and D 1 
jointed together as shown in Fig. 2. The kinetic energy of S' will be the 
same as that in equation (21) less the two terms mentioned above if C, D and E 
are replaced by C', D 1 and E 1 , respectively, and if E' is a fictitious body 
having the mass distribution of C' and D' at the Instant under consideration. 
With perfect rotation of axisymmetrlc rotor R, the inertia dyadic of C and R 
together for Q is fixed C. Therefore a real rigid body C* can have the same 
mass distribution as C and R at all times and have the same motion as C. If 
this choice is made, and in addition, D 1 is chosen to have the mass distribu- 
tion and the motion of D, then E' has the same mass distribution as that of E 
at all times. The kinetic energies K‘ and K , respectively, of S' and the 
original system S consisting of C, D and R can thus be related by 
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K = K' + ^ p^ J q^ + p J q w c ♦ a 


( 22 ) 


It is worth noticing that the common wisdom of simply adding p 2 J to the 
"Inertia of the driven link" to compensate for the drive system dynamics is 
only true when either w c , the angular velocity of the carrier of the drive 
system, or « c • a Is zero. For many robots there are some drive systems that 
do not satisfy either condition. 

The differences between generalized inertia force contributions due to S 
and S' can be worked out based on equation (22). Since generalized inertia 
* 

force F r is related to kinetic energy K by 
c* _ / d 3K 9K v , 

F r " ' ( ~dt W ' r " lf n (2 

T T 

it follows that 


where 


( F r ) s 3 ( F r ) s i + G r r = 1 n 


- p J(pq + a c • a) (q p = q) 


- y 0{q 


9(w c • a) . d 8(« c • a) 9(w c • a). 


+ q C 


dt a^ r 


(q r * q) 


The partial differentiations in equation (25) are most advantageously per- 
formed with C being the reference frame because a is fixed in C. Furthermore, 
If C is the 1th link in an n-link articulated robot as shown in Fig. 3 and D 
is the subsequent link, then it can be shown that 


9(w c • a) { Z| "* a ( ri ^) 

9Q r 0 (r > i) 


9(w C • a) C 3w C > r x z r ) * a (r <; i) 

— - • a = { 

9 ^ r 3q r 0 (r > i) 


i 



(28) 


d? ( z r * *) = ( CwBr x z r ) ' a = ( Bi « Br x z p ) • a 

* 

With equations (26)-(28) In equation (25), one can rewrite G r as 

- yJ (y q. +1 + N a B i -a) (r = i + 1) 

G* = - yJ[(z p • a) q i+1 - q i+1 ( N w Bi x z p ) • a] (r £ i) (29) 

0 (r > i + 1) 


If the generalized inertia forces of system S' have been worked out separa- 
tely, then that of S can be derived using equations (24) and (29). 

As to the generalized active forces, consider that R is acted upon 
through electromagnetic or viscous damping interaction by C and the net result 
of this interaction Is a couple of torque Ta. The laws of dynamics dictate 
that a couple of torque -Ta is also acted on C by R. The contribution of this 
pair of couples to the generalized active forces F r are simply 

yT (r = 1 + 1) 

F_ = (30) 

r 0 (r * 1 + 1) 

No other interaction forces between the bodies in S contribute to the genera- 
lized active forces. For system S', it is easily seen that if a couple of 

torque yTc is assumed to act on D' by C 1 , then the contributions of the 
interaction forces to the generalized active forces are the same as that in 
equation (30). 

Generalization 

In some situations, the joint drive system of a particular joint is 
mounted on the outward link rather than the Inward link of the joint, such as 
that of the second joint of Unimation PUMA robots. The equations derived for 
S can still be applied with D being the inward link and C being the outward 

link. Consider that C is still the ith link, and D is the (i-l)th link. The 

* 

difference terms G r in equation (24) become 


G 


* 

r 


- yJ [y q. + N a Bi 

- yJ [(z f • a) q 1 

0 


a + ^ (z i • a) - ^(V 1 x z.) • a] (r = i) 

q^V 1 x z r ) ♦ a] (r < i) (31) 

(r > i) 


where unit vector a should be in the direction that Is associated in the right 
hand sense with the rotation of R when the joint experiences a positive 
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rotation. The above generalization Is true because all the derivations for 
the generic system In Fig. 1 do not depend on whether C or D is the preceding 
link, with the exception of equations (16) and (19). When D is the link that 
precedes C, unit vectors a and c can be properly changed to maintain the 
validity of these equations. With this in mind, equation (22) remain valid 
for the new case. The difference from the original case Is that w c in 
equation ( 22 ) is a function of q for the present case, and it therefore gives 
rise to the differences between the expressions in equations (29) and (31). 

It is also observed that the above development applies to any axisym- 
metric body designated R that is carried on C and performs fixed axis rotation 
in C. Any gear in a gear train connecting a motor to the link it drives can 
be the rotor R and Its contribution to the change of the generalized inertia 
forces can be identified. It should be noticed that there will be a different 
gear ratio and a different unit vector a for each gear. 

Furthermore, one can see that the only role D plays in equations (29) or 
(31) Is related to the definition of q which is used in equation (16). If the 
generalized coordinates can be properly introduced for the system and the 
angular velocity of R In C can be expressed as a function of these generalized 
coordinates, then the role of D can be eliminated. The above results can thus 
be extended to the drive system for linear joint. They can also be extended 
to complicated gear systems that make up many robot wrist mechanisms such as 
that discussed in [10]. For such a system with three degrees of freedom, 

C R * * 

w = (y 1 q 1 + P 2 q 2 + y 3 q 3 ) a (32) 


where qj, q 2 and q 3 are the generalized coordinates associated with the system 
and ]ii, p 2 and ^3 are the corresponding ratios for R, should be used Instead 
of equation (16). With this, equation (22) is replaced by 

K s = K s + j + y 2 q 2 + P 3 q 3 ) 2 + J^qj + y 2 q 2 + p 3 q 3 ) w c • a 


where S is the system consisting of C and R while S' consists of C' only. 

For a complete n degree of freedom motor-driven robot with speed reduc- 
tions Involved, the above procedure can be applied in the following manner: 

1. Model the system as made up of n rigid bodies connected by the 
appropriate linear or revolute joints with each of these bodies having the 
mass distribution of the actual link and any rotational elements that it 
carries. Formulate the equations of motion for such a model. 

2. For each of the rotational elements, figure out its additional contri- 
butions to the kinetic energy and in turn the contributions to the generalized 
inertia forces as developed above. Add these contributions to the equations 
of motion. 

3. Use equation (30) to work out the generalized active forces. 
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Significance of Reduction Effects 

Dynamic equations for PUMA 560 robot has been explicitly derived with 
parameter values measured or estimated in [11]. Based on the parameter values 
listed, the effects of the drive systems can be estimated. For PUMA 560 
robots, the 2nd and 3rd joint drive systems are mounted on the 2nd link while 
the 4th to 6th drive systems are mounted on the 3rd link. The motors are 
mounted in such a way that their axes of rotation are always perpendicular to 
the 2nd and the 3rd joint axes. If only the inertia matrix, which is the 
congregate of the coefficients of q-; , i= 1 ,.., 6 , and the motor rotor's contri- 
butions are considered, the coefficients to have additional terms include the 
diagonal elements as well as those with indices (1,1), i =2 , . . . , 6, and the off- 
diagonal elements with Indices (4,5), (4,6) and (5,6) due to the coupling of 
the drive systems of the wrist joints. It is understood that the inertia 
matrix is a symmetric matrix, and only the elements in the upper triangular 
part of the matrix are addressed. 

Listed in [11] are inertia contributions of joint drive systems to the 
diagonal elements of the inertia matrix. Here, these are assumed to be mainly 
contributed by the motor rotors In the form of y^j pursuant to equations 
(29) and (30). With this assumption, the constant coefficients of the domi- 
nant terms, involving sine and cosine functions of joint angles, of the matrix 
elements effected can be compared with those additional contributions propor- 
tional to yJ as computed based on equations (29) and (30). Table 1 shows the 
list. 


Table 1 Effects of Drive System on Inertia Matrix 


Element 

1,1 

2,2 

3,3 4,4 

5,5 

6,6 

1,2 1,3 

1,4 

1,5 

1,6 

Dorn. Coef. 

2.57 

6.79 

1.16 .20 

.18 

.19 

.69 .13 

1.64 

1.25 

4. 

as in [2] 







E-3 

E-3 

E-5 

Coef. of 

1.14 

4.71 

.83 .2 

.18 

.19 

.04 .015 

2.60 

2.5 

2.5 


add. term E-3 E-3 E-3 

2 

It can be seen from Table 1 that y J terms are dominant In the diagonal 
elements of the inertia matrix. They are included in the equations in [11] 
while the additional contributions to the off-diagonal elements that are 
proportional to yJ are not included. Due to the fact that their contributions 
remain constant when the robot posture is changed, the percentage variations of 
the diagonal elements is smaller than what it would be if the robot is a direct 
drive one. This makes fixed gain control more likely to succeed. Although some 
of the latter ones are dominant, they are still very small compared to the 
(1,1), or even the (4,4), (5,5) or (6,6) element. As to the elements (4,5), 
(4,6) and (5,6), since the coupling relationship between the drive systems are 
not discussed in [11], the additional contributions cannot be estimated. It is 
reasonable to predict that they will be more significant than those for other 
off-diagonal elements judging from equation (33). 

Contributions due to the motor rotors to those terms second order in qj, 
1=1,.., 6, can also be estimated. Again, some of them may be dominant in 
the coefficient of a particular term, but their effect to the complete 
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equations may not be significant unless qj, 1=1,.., 6, assume significantly 
high values. 

Conclusions 

Although it has been known to joint drive system designers that the 
Inertia properties of the motor rotor and other elements connected to it are 
important factors in determining the system dynamic response, it has not been 
elaborated in so many articles on robot dynamics. The contributions in the 
form of ysJ in the inertia matrix can dominate some of the matrix elements. 

Other contributions proportional to pJ are less significant, but they may not 
be negligible in all cases. With more and very different robots to be developed 
in the future, it is important to know what need to be included in the dynamic 
model for it to have sufficient fidelity. The methodology set forth in this 
paper provides a means to gain the necessary information for a sound judgment. 
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